Exact dynamics of entanglement and entropy in structured environments 
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We study the exact entanglement dynamics of two qubits interacting with a common zero- 
temperature non-Markovian reservoir. We consider the two qubits initially prepared in Bell-like 
states or extended Werner-like states. We study the dependence of the entanglement dynamics on 
both the degree of purity and the amount of entanglement of the initial state. We also explore the 
relation between the entanglement and the von Neumann entropy dynamics and find that these two 
quantities are correlated for initial Bell-like states. 
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I. INTRODUCTION 



Entanglement is a key resource in modern quantum in- 
formation theory and technology. Entangled states play a 
central role in quantum key distribution, superdense cod- 
ing, quantum teleportation and quantum error correction 
However, realistic quantum systems are never com- 
pletely isolated from their surroundings. The inevitable 
interaction between a system and its environment leads 
to decoherence phenomena and degradation of entangle- 
ment 0. 

Once entanglement has been lost, it cannot be restored 
by local operations 0]. It is therefore important to un- 
derstand the process of disentanglement in order to con- 
trol the effects of noise and to preserve entanglement. 
For these reasons the study of entanglement dynamics 
of quantum systems in realistic situations has become of 
increasing importance. 

Recently much attention has been devoted to the pro- 
cess of finite-time disentanglement, also known as "en- 
tanglement sudden death" (ESD), in a bipartite system 
This phenomenon consists in the complete disap- 
pearance of the bipartite entanglement in a finite time, 
as opposed to the time-evolution of the coherences of the 
single parts which vanish only asymptotically. 

ESD sets a limit on the life-time and usability of entan- 
glement for practical purposes. Hence, lots of efforts have 
been done in order to understand the conditions under 
which ESD occurs Q. Finite-time disentanglement has 
been found in different physical systems, such as qubits 
or harmonic oscillators [6|. ESD appears when the two 
parts of the system interact with either independent en- 
vironments [1, 0] or a common one @ . The first studies 
within the Markov approximation have been extended to 
non-Markovian environments, where the memory of the 
reservoir adds revivals to entanglement dynamics [To| . 

In Ref. [9( we have studied entanglement sudden death 
when two-qubits are prepared in a Bell-like state with two 
excitations, and entanglement sudden birth for qubits 
prepared in a separable state. Here we focus on the 
dynamics of a class of states having an "X" -structure 
density matrix, namely the extended Werner-like states 
(EWL). This class of states plays a crucial role in many 
applications of quantum information theory, such as tele- 



portation and quantum key distribution [12(]. More- 
over, such a choice will give us the chance to study how 
the entanglement dynamics and its revivals are related to 
the purity and the amount of entanglement of the initial 
state. We also focus on the interplay between entangle- 
ment and mixedness for the qubits states. We investi- 
gate the connection between these two quantities, com- 
paring concurrence and von Neumann entropy dynamics 
for initial Bell-like states, and finding clear correlations 
between them. 

The paper is organized as follows. In Sec. II we re- 
view the exactly solvable model of two qubits interact- 
ing with a common Lorentzian structured reservoir. In 
Sec. Ill we study the entanglement and von Neumann 
entropy dynamics for initial Bell-like states and we prove 
that sudden death of entanglement can never occur if the 
qubits are initially in a mixed state having at most one 
excitation. In Sec. IV we focus on the time evolution 
of EWL states, and compare our results with those ob- 
tained in Ref. [13j | for independent structured reservoirs. 
Finally, we summarize our results in Sec. V. 



II. THE MODEL 

In this section we describe the model we use to study 
the dynamics of two two-level systems (qubits) interact- 
ing with a common zero-temperature bosonic reservoir. 
Our approach is non-Markovian and non-perturbative, 
i.e., it does not rely on either the Born or the Markov 
approximations @. 

The Hamiltonian of the system, in the rotating wave 
approximation, is given by H = Hq + Hi nt , 



flint = 0+ + cr 



k 



h.c. 



(1) 



(2) 



where a± and a± are the Pauli raising and lowering op- 
erators for qubit A and B respectively, luq is the Bohr 
frequency of the two identical qubits, a k and at , u> k and 



2 



<7fc are the annihilation and creation operators, the fre- 
quency and the coupling constant of the field mode k, 
respectively. 

In order to solve the dynamics of the two qubits we 
need to specify the properties of the environment. In 
the following we assume that the two qubits interact 
resonantly with a non-Markovian Lorentzian structured 
reservoir, such as the electromagnetic field inside a lossy 
cavity [14| , having spectral distribution 



J(u) 



n 2 



2ir (w - lo ) 2 + (r/2) 2 



(3) 



where T is the width of the Lorentzian function and f2 
the coupling strength. 

We are mainly interested in the dynamics of entangle- 
ment and in the effects that the non-Markovian reser- 
voir induces on the correlation between the two qubits. 
To quantify entanglement we use the Wootters concur- 
rence [TB|, defined as C(t) = max{0, \Ai — — — 
V^A^}, where {A^} are the eigenvalues of the matrix 



R 



Op- 



{<Jy ®cry), with p* denoting the com- 

A/B 



plex conjugate of p and o y are the Pauli matrices for 
atoms A and B. This quantity attains its maximum value 
of 1 for maximally entangled states and vanishes for sep- 
arable states. 

We focus now on the dynamics of initial "X" states. 
We use the method described in Ref. § to calculate the 
time evolution of the density matrix 



P(t) 




(4) 



which is written in the basis {|00), |10), |01>, |11)}. Due 
to the structure of the differential equations for the den- 
sity matrix elements [see Eqs. (5) and (6) in Ref. @], the 
"X" form is preserved during the evolution. In the Ap- 
pendix we present the analytical solution in the Laplace 
transform space for a particular type of "X" state. 

For this class of states the concurrence assumes a sim- 
ple analytic expression 



where 



C(t) =max{0,C 1 (t),C 2 (t)}, 



Ci(t) = 2\w(t)\ - 2y/b(fr(fi, 
C 2 (t) = 2\z(t)\-2y/a(t)d(t). 



(5) 



(6) 



We notice that coherences give a positive contribution to 
Ci(t) and C2(t) and so to concurrence, while the negative 
parts involve populations only. 

In the next section we will also consider the evolution of 
the mixedness of the two qubit state, which we quantify 
through the von Neumann entropy, defined as 



Von Neumann entropy is equal to zero for pure states, 
and attains its maximum value (equal to In N with N the 
dimension of the Hilbcrt space) for a maximally mixed 
state. 



III. ENTANGLEMENT AND VON NEUMANN 
ENTROPY DYNAMICS FOR BELL-LIKE STATES 

Here, we seek for the connection between the dynamics 
of entanglement and von Neumann entropy of two qubits 
prepared in Bell-like states 



and 



|$) =o|10) +e M (l -a 2 ) 1/2 |01) 



I*) =a|OO)+e 40 (l-a 2 ) 1 / 2 |ll). 



(8) 



(9) 



S(p) = -Tr{p(t)Hp(t))}. 



(7) 



Our aim is to understand the interplay between these 
two different physical quantities, in particular when pe- 
culiar phenomena such as ESD or "entanglement sudden 
birth" (ESB) dd QjJ , and revivals of entanglement oc- 
cur. 

The entanglement dynamics of two qubits in a 
Lorentzian structured reservoir, prepared in a Bell-like 
state with two excitations as in Eq. ©, has been pre- 
sented in Ref. [9J. There we have studied in detail the 
difference with the common Markovian reservoir case and 
the independent reservoirs non-Markovian case. The evo- 
lution of entanglement for a Bell-like state as in Eq. ([5]) 
has been studied also in Ref. [18|. 

For both the Bell-like states in Eqs. © and © the 
entanglement dynamics is the result of two combined ef- 
fects: the backaction of the non-Markovian reservoir and 
the reservoir-mediated interaction between the qubits. 

The memory effects due to the non-Markovianity of 
the reservoir causes oscillations in entanglement dynam- 
ics. These oscillations are typical also of the independent 
reservoirs case [T(|, but they disappear completely for 
Markovian reservoirs Q. The sharing of the reservoir 
plays also a special role. Indeed, the common reservoir 
provides an effective coupling between the qubits, and 
so consequently creates quantum correlations between 
them. As a result, qubits prepared in a factorized state 
can become entangled due to the interaction with the 
common reservoir. This is in contrast with the indepen- 
dent reservoirs case in which a factorized state of the 
qubits can never evolve into an entangled state. 

The results in Ref. [l8| show that ESD does not oc- 
cur for a Bell-like state with one excitation as in Eq. ([8]) 
for any value of a 2 . ESD does not appear for the same 
Bell-like state even if the dipolar interaction between the 
qubits is included [T9j |. Actually, a straightforward cal- 
culation shows that for every pure or mixed state of the 
qubits containing at most one excitation, ESD and ESB 
cannot take place. In fact, the density matrix describ- 
ing a generic mixed state with maximum one excitation, 
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written in the same basis of Eq. (|4|) , has the form 



P(t) = 




(10) 



The expression for the concurrence for any value of the 
parameters is 



C(t) = max{0,2|z(t)|}. 



(11) 



Here concurrence is directly given by the coherence be- 
tween the 1 10) and |01) states. Since the coherence van- 
ishes in asymptotic way, there cannot be ESD for any 
generic state with maximum one excitation. Analogously, 
entanglement can be smoothly generated but it cannot 
suddenly appear. This result does not depend on the de- 
gree of purity of the state. This is true as long as the form 
of the density matrix in Eq. (flO)) is maintained. On the 
other hand, if some population is transferred to the two 
excitations state then ESD can appear. This is the case 
of two qubits in a Bell state interacting with a non-RWA 
common reservoir [20| . 

One problem we want to address here is to understand 
how the amount of entanglement changes with the pu- 
rity of the state. In the next section we will explore this 
aspect by preparing the qubits in a particular class of 
mixed states. Furthermore, it is useful to see how the 
mixedness of the initial state changes with time because 
of the interaction with the environment, and if there ex- 
ists some connection between the von Neumann entropy 
and entanglement dynamics. 

In order to provide a clear interpretation of the dy- 
namics, it is useful to recall the four-state effective de- 
scription Indeed, the state of the two qubits is 
equivalent to a four-state system, in which three states 
{|00), |+) = (| 10) + |01»/V2, |H)} interact in lad- 
der configuration with the electromagnetic field, and the 
fourth state {|-) = (|10) - |01))/V2} is completely de- 
coupled from the other states and the electromagnetic 
field. The states |+) and |— ) are known as super-radiant 
and sub-radiant states, respectively. 

We look at the evolution of the von Neumann entropy 
of the qubit pair when the state is initially prepared in a 
Bell-like state of the form (JSJ and . It is particularly 
interesting to see how the degree of purity of the state 
evolves when the qubits undergo ESD. Since ESD occurs 
when a 2 < 1/4 |{|, we choose a 2 — 1/20 and 8 — for 
the sake of convenience. In Fig. [1] we notice that ESD ap- 
pears for high value of mixedness of the system (entropy 
peaks coincide with minima of concurrence). Moreover, 
the revivals of entanglement appear roughly in correspon- 
dence of the minimum of the von Neumann entropy, when 
the state becomes purer. A careful analysis shows that 
the dynamics of entropy follows the population of the 
super-radiant state p++(t), (in the Appendix we provide 
the analytic expression of the von Neumann entropy for 
this initial state). As Fig. [T] shows, they attain their 
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FIG. 1: (Color online) Dynamics in a common Lorentzian 
structured reservoir as a function of scaled time for two atoms 
prepared in the Bell- like state <& with a 2 = 1/20 and 6 = 0. 
Solid red line is concurrence; dotted-dashed blue line is von 
Neumann entropy; dotted green line is the population of the 
super-radiant state p++(t). In the inset: dotted green line is 
the super-radiant state population; dotted-dashed black line 
is the ground state population; solid light blue line is the 
excited state |11) population. 




FIG. 2: (Color online) Dynamics in a common Lorentzian 
structured reservoir as a function of scaled time for two atoms 
prepared in the Bell-like state $ with a 2 = 1/2 and 6 — 0. 
Solid red line is concurrence and dotted-dashed blue line is 
von Neumann entropy, whiich overlaps with the super-radiant 
state population. 



relative maxima and minima at the same times. More- 
over, we can deduce from Eqs. ([5]) and ^ that entan- 
glement dynamics is a function of the population of the 
super-radiant state, which is p++(t) = 2y / b(t)c(t) for the 
particular initial state of Eq. ([9]). Therefore whenever 
the population p ++ {t) reaches its relative maxima, the 
state attains a maximum value of mixedness, the time- 
dependent part of the concurrence (2|w(t)| — p ++ (t)) be- 
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FIG. 3: (Color online) Dynamics in a common Lorentzian 
structured reservoir as a function of scaled time for two atoms 
prepared in the Bell-like state $ with a 2 = 1/5 and 6 — 0. 
Solid red line is concurrence; dotted-dashed blue line is von 
Neumann entropy; dotted green line is super-radiant state 
population. In the inset: dotted green line is super-radiant 
state population; dotted-dashed black line is ground state 
population; dashed pink line is sub-radiant state population; 
solid yellow line is absolute value of the coherence between 
super-radiant and sub-radiant states. 




FIG. 4: (Color online) Dynamics in a common Lorentzian 
structured reservoir as a function of scaled time for two atoms 
prepared in the mixed state (j 12|) with a 2 — 0.75. Solid red 
line is concurrence; dotted-dashed blue line is von Neumann 
entropy; dotted green line is super-radiant state population. 

comes negative and entanglement disappears. On the 
other hand, whenever the population of the super-radiant 
state reaches a minimum, the population of the |11) ex- 
cited state and the |00) ground state have their maxima 
(see inset in Fig. [T|), and the system goes toward a Bell- 
like state \^f). As a consequence the system becomes 
purer and entanglement is partially recovered. 

When the system is prepared in a Bell-like state as the 



one of Eq. ([8]) there is no ESD, however it is still of inter- 
est to understand how the revivals of entanglement are 
related to the degree of purity of the state. For the Bell 
state |+) entanglement has exactly the same dynamics 
of the population of the super-radiant state, as shown in 
Fig. [2 Moreover the zeroes of entanglement and entropy 
coincide. For those times, in fact, the system goes into 
the ground state which is pure and factorized. When 
some population returns in the super-radiant state en- 
tanglement is recovered, and the state is again mixed. 

When a 2 ^ 1/2 and/or 8^0 the initial state is a su- 
perposition of super-radiant and sub-radiant states. Al- 
though the sub-radiant state does not evolve in time, be- 
ing decoupled from the super-radiant and ground states, 
it affects the entanglement and entropy dynamics. For 
the sake of convenience we choose a 2 = 1/5 and 9 = 0. 
In Fig. [3] we see that entropy still follows closely the time 
evolution of the super-radiant state population, having 
its relative minima in the same positions of the zeroes 
of p ++ {t), (the analytic expression of the von Neumann 
entropy is in the Appendix). On the contrary, entan- 
glement has new relative maxima when the population 
of the super-radiant state is zero. This is due to the 
presence of the sub-radiant state. In fact, in this case 
both the super-radiant and the sub-radiant states con- 
tribute to the total entanglement. Thus there are two 
different sets of entanglement maxima, those associated 
with the maxima of the super-radiant state population, 
and those associated to the sub-radiant state. Entan- 
glement is zero whenever the population of the super- 
radiant state p++{t), the population of the sub-radiant 
state p (t) , and the absolute value of the coherence be- 
tween super-radiant and sub-radiant states /?_| (t), are 

equal. This can be explained in the light of the expres- 
sions in Eqs. (|A.1[) in the Appendix, where b(t), c(t) 
and z(t) are written as a function of those quantities. 

In fact, when p ++ (i), p (t) and \p^ (t)\ have equal 

value, then z(t) is equal to zero, irrespective of the sign 

of p^ (t). Specifically, this happens in correspondence 

to the 1st, 4th and 5th zeroes of entanglement, where 
b(t) = and c(t) = 2p__ = k. In this case the state 
of the system becomes (1 - fc)|00)(00| + fe|01)(01|, which 
is clearly not entangled. Similarly, for the 2nd and 3rd 
zeroes, b(t) = k, c(t) = and the non-entangled state is 
(1- fe)|Q0> <00] +fc|10)(10|. 

To conclude this section we consider the following 
mixed and factorized initial state 

p(0) = (a 2 \0 A )(0 A \ + (l-a 2 )\l A )(l A \) 

®(a 2 \0 B )(0 B \ + (l-a 2 )\l B )(l B \). 

The dynamics of this state in a common non-Markovian 
reservoir is characterized by ESB and revivals of disen- 
tanglement, as we have shown in Ref. 9]. In Fig. Q] we see 
how these interesting features are related to the degree of 
purity of the state. As we have seen before, the positions 
of maxima and minima of the entropy and of the super- 
radiant state population match. The sudden creation 
of entanglement happens roughly when the entropy hits 
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its first minimum. Entanglement is again lost when the 
amount of mixedness increases, and it reappears again 
when the entropy reaches another minimum. 

IV. WERNER STATE 

In this section we study the two qubits entanglement 
dynamics of extended Werner-like states in a common 
zero-temperature Lorentzian structured reservoir. Our 
goal here is to study how, starting from an initial state 
that is not perfectly pure, as in realistic experimental 
conditions, affects the entanglement dynamics, and in 
particular the occurrence of ESD and ESB phenomena. 




A. Terminology and previous works 

The standard two-qubit Werner state, introduced in 
1989 by Werner [2l|, is defined as 

Pw=r\-)(-\ + ^l, (13) 

where |— ) is the singlet state. In Ref. [21] Werner demon- 
strated that while pure entangled states always violate 
the Bell inequality, mixed entangled states might not. 
The Werner state is the first entangled state to be proven 
not violating any Bell inequalities [2l|. The generalized 
or Werner-like states are defined as 

p W L = r\M){M\ + ^^I (14) 

with \M) one of the four maximally entangled Bell states. 
For a given r, Werner and Werner-like states exhibit 
the same entanglement. A further generalization is the 
extended Werner- like states (EWL), containing a non- 
maximally entangled state part, which are defined as 

P%WL = r\$)($\ + ^I, (15) 

where |$) = a|10) + e w {1 - a 2 ) 1/2 |01), and 

pt WL =r\*){*\ + ^I, (16) 

with |*> = a|G0) + e l6 {l - a 2 ) 1/2 |ll). 

The dynamics of entanglement of two qubits prepared 
in Werner, Werner-like or extended Werner-like states 
has attracted a lot of attention. The appearance of 
ESD has been demonstrated for two qubits prepared in 
a Werner-like state and interacting with Markovian in- 
dependent reservoirs [22], [23j or with independent noisy 
channels [24[. Violation of Bell inequality has been ex- 
amined in independent thermal reservoirs [25l |2(| . En- 
tanglement has been studied in a Markovian common 
thermal reservoir [27|], in the presence of collective de- 
phasing [28| and in the maximal noise limit [29[ . Finally, 



FIG. 5: (a) Concurrence as a function of scaled time and r 
for two atoms prepared in the Werner state (I15p ($) with 
a 2 = 1/2 and 9 = 0. 

the dynamics of an EWL in two independent Lorentzian 
structured reservoirs has been investigated in Ref. 

Moreover, Werner and Werner-like states have been 
used so far in many applications in quantum informa- 
tion processes such as teleportation [llj and entangle- 
ment teleportation [30]. The experimental preparation 
and characterization of the Werner states have also been 
widely investigated. Werner states are prepared via spon- 
taneous parametric down-conversion [3l| or using a uni- 
versal source of entangleme nt 32l| . and used in ancilla- 
assisted process tomography |33j and secure quantum key 
distribution [T3 |. 

B. Entanglement dynamics 

We evaluate the dynamics of the entanglement of EWL 
states as a function of the initial amount of mixedness, 
controlled by the purity parameter r, and as a function 
of the initial degree of entanglement measured by a 2 . 

Looking at the evolution of this kind of states gives 
the possibility to study how the degree of purity of the 
initial state influences the entanglement dynamics. We 
also compare our results with those obtained in Ref. [l3[ 
for two qubits in two independent Lorentzian structured 
reservoirs. 

Figures [51 [5] and [7] show the entanglement dynamics 
as a function of the dimensionless quantity jot (with 
7o = 4£! 2 /r the Markovian decay rate of the atoms) 
and of the purity parameter r, for the two EWL states 
in Eqs. (fl~o| and (flUj) . In Fig. [5j the qubits are ini- 
tially in the state given by Eq. (fTo| with a 2 = 1/2 and 
6 = 0. When r = 1 the qubits are prepared in the super- 
radiant state, entanglement exhibits oscillations and ESD 
is never present. However, whenever a little amount of 
mixdness is added, periods of finite-time disentanglement 
immediately occur. This is due to the appearance of some 
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FIG. 6: (a) Concurrence as a function of scaled time and a 2 
for two atoms prepared in the Werner state (|15p (<3>) with 
r = 1/2 and 9 = tt. 




FIG. 7: (a) Concurrence as a function of scaled time and r 
for two atoms prepared in the Werner state (I16|l with 
a 2 = 1/2 and = 0. 

population in the excited state 1 11). As a consequence the 
oscillating part of the concurrence {2|z(t)| — 2y/a(t)d(i)} 
can become negative. Due to the non-Markovianity of the 
system and to the effective coupling provided by the com- 
mon reservoir, ESD regions are followed by revivals, and 
for long times entanglement reaches a stationary value. 
When r < 1/3 the state is initially factorized, but as 
time passes, because of the reservoir-mediated interac- 
tion, entanglement between the qubits is suddenly cre- 
ated. The non-Markovianity of the reservoir enriches the 
dynamics causing eventually revivals of disentanglement. 
Eventually the entanglement reaches a non-zero station- 
ary value. The amount of entanglement that has been 
created depends on the population of the sub-radiant 
state, which carries the entanglement. 

If the qubits are prepared in the state given by Eq. JT5|) 



with a 2 — 1/2 and 9 — tt, the dynamics is completely dif- 
ferent compared to the previous 9 = shown in 
Fig. [SJ In fact when r = 1 the qubits are prepared in the 
sub-radiant state and concurrence does not evolve, being 
at any time equal to 1. When the mixed part is present 
little oscillations appear, but eventually the concurrence 
attains the stationary value (l+3r)/4. When r < 1/3 the 
initial state is factorized, entanglement is suddenly cre- 
ated, a revival of disentanglement is present, and again 
after some oscillations concurrence reaches its stationary 
value. 

When a 2 ^ 1/2 and/or 9 ^ 0, tt the non- mixed part of 
the initial state is a superposition of super-radiant and 
sub-radiant states. Thus the entanglement dynamics de- 
pends on the weights of those two states. As a result 
there is a wide variety of entanglement dynamics in be- 
tween the two asymptotic behaviors described above. 

In Fig. [7] the qubits are prepared in the state given by 
Eq. (JTBj! with a 2 = 1/2 and 9 = 0. Note however that 
anyway the results are independent of the choice of the 
relative phase 9. For r = 1 the initial state is the Bell 
state ( 1 00) 4- |ll))/\/2. For this initial condition not only 
there is no ESD, but also concurrence vanishes only for 
infinite time. When an increasing amount of mixedness 
is present in the initial state, finite-time disentanglement 
appear. ESD is then followed by revivals and, as ex- 
pected, a certain amount of entanglement is preserved. 
Analogously to the other EWL state, for r < 1/3 the ini- 
tial state is factorized. For the same reasons previously 
mentioned, entanglement is suddenly created, momentar- 
ily deteriorated, and it finally goes to the stationary value 
r/4, coinciding with the population of the sub-radiant 
state. For this type of EWL state different choices of a 
lead to the same qualitative behavior of the entanglement 
dynamics. 

The results we just described for two qubits in a com- 
mon structured reservoir are quite different from those 
presented in Ref. 13] for two qubits in two indepen- 
dent Lorentzian structured reservoirs. First of all, if two 
uncoupled qubits interact with two independent reser- 
voirs entanglement cannot be created from a factorized 
state. Hence, the ESB region, characterizing the dynam- 
ics in a common reservoir, is not present in the results 
in Ref. [l3|], for any of the EWL states. Moreover, for 
qubits prepared in the initial state of Eq. (|16p . we no- 
tice that the reservoir-mediated interaction between the 
qubits keeps the value of concurrence higher compared 
to the two independent reservoirs case. 

The crucial difference between the common and inde- 
pendent reservoirs cases is that for qubits in two inde- 
pendent reservoirs the decoherence-free sub-radiant state 
does not exist. When the qubits are in two independent 
reservoirs, the entanglement dynamics of the |+) and |— ) 
states are the same. As a consequence the relative phase 
6 in Eqs. ([8]) and (fT5|) does not affect the results. This 
is definitely in contrast with the results we present in 
Figs. O and [51 showing two completely different asymp- 
totic behavior when changing the relative phase 9 of a 7r 
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factor. 



V. CONCLUSIVE REMARKS 



stand the interplay between entanglement and entropy, 
and that the present results shed new light on their dy- 
namical relation. 



In this paper we have investigated the connection be- 
tween entanglement and entropy dynamics in a system 
of two qubits interacting with a common zero tempera- 
ture non-Markovian reservoir. We have used the exactly 
solvable model presented in Ref. [9j , in the Appendix we 
attach the analytical solution in the Laplace transform 
space in the case of EWL states. 

We have compared the entanglement and von- 
Neumann entropy time-evolution for two qubits prepared 
in a Bell-like state. We have noticed that ESD seems to 
appear when the state becomes highly mixed, whereas re- 
vivals of entanglement are associated to minima of the en- 
tropy, where the state becomes purer. On the other hand, 
when starting from a factorized state, sudden birth of 
entanglement occurs for lower values of mixedness, while 
revivals of disentanglement are accompanied by peaks of 
the entropy. 

For a Bell-like state with one excitation the picture 
seems to be more complex. However, by studying the 
dynamics of the population of the super-radiant and sub- 
radiant states, we have realized that two sets of maxima 
of entanglement can be identified. The super-radiant 
state mainly carries the entanglement when it is max- 
imally populated; when its population is zero, the en- 
tanglement is associated to the sub-radiant state. When 
the populations of the states are equal, the entanglement 
vanishes. 

We have also considered the entanglement dynamics of 
a particular class of mixed states, the extended Werner- 
like states. We have demonstrated that the amount of 
purity of the initial state plays a key role in the entangle- 
ment dynamics, controlling the appearance of ESD and 
ESB phenomena. For stronger non-Markovian conditions 
the dynamics exhibits stronger and longer lasting entan- 
glement oscillations and an increasing number of dark 
periods and revivals as well. On the other hand, in the 
Markovian regime, no oscillations are present, however 
the basic features of the dynamics, and so the ESD and 
ESB regions, are still present. 

In general, entanglement and mixedness are two differ- 
ent physical quantities characterizing the degree of non- 
classicality of a quantum state. We think that it is im- 
portant also from fundamental point of view to under- 
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APPENDIX 

Here we present the exact analytic solution for two 
qubits interacting with a Lorentzian structured reservoir, 
when the qubits are prepared in an EWL state as in 
Eqs. (fT5|) and (TIT))) . We provide the expressions for the 
density matrix element in Eq. Q as a function of the 
solution pij of the pseudomode master equation [34| in 
Eqs. (5) and (6) of Ref. [| 



a(t) 


= Paa(t) + 


Pbb(t) + 


Pcc(t), 






b(t) 


_ P++(*) - 


f p— (t) 


+ P+- 


(t) + p_. 








2 






c(t) 


_ P++(*)- 


fp— (t) 


-P+- 


W-p- 








2 






z(t) 


_ P++(*)- 


~P-(t) 


- P+- 


W + P- 


f(t) 






2 






d(t) 


= 1 - a(t) 


-&(*)- 


c{t) = 


p//(*)> 




w(t) 


= Paf(t), 










p++(t) 


= Pdd(t) + 


Pee(t). 









(A.l) 



The subscripts a, b and c refer to the states in which 
both the qubits are in the ground state, and the pseudo- 
mode has 0, 1 or 2 excitations respectively. In d and e 
the atomic system is in the super-radiant state, and the 
pseudomode has or 1 excitations. / is the state having 
excitations in the pseudomode and both the qubits in 
their excited states. 

The analytic expressions are given in the Laplace 
transform space. In the following we provide the solu- 
tions for a EWL state of the form ([TBI : 



J 



Paa{t)=C- 1 [-^- s {-l+r + 



64r 2 ft 4 (-l + r)l(s) 8m 2 (l + r + 4a Vl - a 2 rcos( 



j(s)k(s) 



(A.2) 



i r 2ft 2 , 



8rn 2 {-i + r)i( s ) 



4a\/r 



err cos 



»}■ 



(A.3) 
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p cc (i)=£ { ^ }, (A.4) 

Pddit) = c- l i——{-—^{srn 2 (-i + r)(6r 5 + 3ir 4 s + 2or 3 (-2ft 2 + 3s 2 ) + sr 2 s(-20ft 2 + lis 2 ) 

L4j(s) v fe(s) 

+ 8r(56fi 4 - — fl 2 s 2 + 3s 4 ) + 4(120ft 4 s - 2ft 2 s 3 + s 5 ))) + (8ft 2 + (r + s)(T + 2s))(l + r + 4a a/I - a 2 r cos0))}, 
2 (A.5) 



,._ t r 2ft 2 (-l + r)(6r 3 + 8m 2 + 13r 2 s + 12ft 2 s + 9rs 2 + 2s 3 )~i 

( ^ ), (A.6) 



p ff (t) ^^j-—— ((-1 + r )( 6r5 + 3ir4s + 4r 3 (38ft 2 + 15s 2 ) + T 2 s(412ft 2 + 55s 2 )+ 



4fc(s) 

8r(40ft 4 + — fi 2 s 2 + 3s 4 ) + 4(72ft 4 s + 26ft 2 s 3 + s J 



(A.7) 



x f r(T + 2s)(-l + 2a 2 + 2WT^? S infl)) ^ 
I 2(4ft2 + s(r + 2s) /' (A - 8) 

I 

p_+(t) = pi_(t), (A.9) j(s) = (r + 2s)(8r! 2 + s(r + s)), (A.13) 



rfl — 2a\/l — a 2 cos 9) 1 — r , k . 
P—(t) = - \ (A ' 10) 

Paf(t) = 0. (A. 11) 

tv. w f k + wv, t7wt ++ v. Z(s) = (6r 3 + 3ir 2 s + T(56ft 2 + 45s 2 ) + 20(3ft 2 s + s 3 )). 

lhe solutions for both the hyvL states are written as f \ \a\ 

functions of k(s), j(s), l(s) defined as 

k(s) — (16rft 2 + 2r 2 s + 24£l 2 s + 3I\s 2 + s 3 ) Here we present the solutions for a EWL state of the 

( 3 r 3 + 28 m 2 + nr 2 s + 24ft 2 s + 121V + 4s 3 ), (A ' 12) form m : 



m r-^- l ( io. i 2 , 8rr> 2 (~l + r) 64r 2 ft 4 (-l + (-3 + 4^p)r)/(s) U 

**> = ^{^(1 -r- gfe^M )}, (A, 6 ) 
M«)=^{- ^- 1 + ^ + 4Qa ) r ^ + a ) }, (A.17) 



3ir 4 s 



M*) = -C" 1 {- (( 8fi2 + ( r + s X r + 2s ))(- 1 + r ) + fcTJ) (8rf! 2 (-l + (-3 + 4a 2 )r)(6r 5 + 
+20r 3 (-2ft 2 + 3s 2 ) + 5r 2 s(-20ft 2 + lis 2 ) + 8r(56ft 4 - y^ 2 « 2 + 3s 4 ) + 4(120ft 4 s - 2ft 2 s 3 + s 5 )))) j 



(A.18) 
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Pee(t)=C ( ^ }, (A.19) 

Pf f (t) =C- l {--^-—{{-\ + (-3 + 4a 2 )r)(6r 5 + 3ir 4 s + 4T 3 {38fl 2 + 15s 2 ) + r 2 s(412f7 2 + 55s 2 )+ 

% . (A.20) 

8r(40fi 4 + yfi 2 * 2 + 3s 4 ) + 4(72ft 4 s + 26Q 2 s 3 + s 5 ))) }, 

- m _ r -i r g^VlEgP! + ( r + *)(r + 2 *)) i rA91l 

pa fW- L l 4ra 2 + r 2 s + i2fi 2 + 3rs 2 + 2s 3 /' (A ^ ij 

i 

p (t) = p+_(t) = £-+(*) = 0. (A.22) entropy for the Bell-like states in Eqs. © and © 

We also provide the expressions of the von Neumann 
I 



=~ («*) + d(*))£o 5 [4] - 2/5 ++ (<)L .g[ P++ (t)] - (o(t) + d(t) - f(t))(Log[a(t) + d(t) - /(*)]) 
^ (A. 23) 

- (o(t) + d(t) + /(i))Lop[a(t) + d(t) + /(t)]) , 

^ $ =\(-2a{t)Log[a(t)] - (p ++ (t) + p__(t) - j(t))Log[p ++ (t) + /3__(t) - 

2 v (A.24) 

- (p++(t) +p__(t)+j(t))io fl [p ++ (t) +p__(t)+j(*)])- 
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